function exp_log = dirichletln_exp(a)

exp_log = digamma(a) - digamma(sum(a));
